Pion propagation in the linear sigma model at finite temperature 



Alejandro Ayala, Sarira Sahu 
Instituto de Ciencias Nucleares 
Universidad Nacional Autonoma de Mexico 
Aptartado Postal 70-543, Mexico Distrito Federal 04510, Mexico. 

We construct effective one-loop vertices and propagators in the linear sigma model at 
finite temperature, satisfying the chiral Ward identities and thus respecting chiral symmetry, 
treating the pion momentum, pion mass and temperature as small compared to the sigma 
mass. We use these objects to compute the two-loop pion self-energy. We find that the 
perturbative behavior of physical quantities, such as the temperature dependence of the 
pion mass, is well defined in this kinematical regime in terms of the parameter m^/iiT 2 f 2 
and show that an expansion in terms of this reproduces the dispersion curve obtained by 
means of chiral perturbation theory at leading order. The temperature dependence of the 
pion mass is such that the first and second order corrections in the above parameter have 
the same sign. We also study pion damping both in the elastic and inelastic channels to this 
order and compute the mean free path and mean collision time for a pion traveling in the 
medium before forming a sigma resonance and find a very good agreement with the result 
from chiral perturbation theory when using a value for the sigma mass of 600 MeV. 

PACS numbers: ll.10.Wx, 11.30. Rd, 11.55.Fv, 25.75.-q 



I. INTRODUCTION 

A basic ingredient for the understanding of several physical processes taking place in hadronic plasmas are 
the propagation properties of pions. The relevant scenarios include dense astrophysical objects such as the 
cores of neutron stars and the evolution of the highly interacting region formed in the midsts of high-energy 
heavy-ion collisions. In the latter context, pions are the most copiously produced particles and their (mainly 
attractive) interactions within a dense nuclear environment have been speculated to give rise to interesting 
collective surface phenomena Q. Moreover, it is well known that the pion group velocity is an important 
piece of information necessary to properly account for the component of the dilepton spectrum originated 
in the hadronic phase, below the critical temperatures and densities for the formation of a quark-gluon 
plasma 0. 

The hadronic degrees of freedom are usually accounted for by means of effective chiral theories whose 
basic ingredient is the fact that pions are Goldston bosons coming from the spontaneous breakdown of chiral 
symmetry. Chiral perturbation theory (ChPT) is one of such effective theories that has been employed to 
show the well known result that at leading perturbative order and at low momentum, the modification of 
the pion dispersion curve in a pion medium is just a constant, temperature dependent, increase of the pion 
mass ||. In spite of the mounting complexity introduced by the number of couplings required at next to 
leading order, ChPT has also been used in a two-loop computation of the pion self-energy Q and decay 
constant ||. A striking result obtained from such computations is that at second order, the shift in the 
temperature dependence of the pion mass is opposite in sign and about three times larger in magnitude than 



the first order shift, already at temperatures on the order of 150 MeV. This result might signal either the 
breakdown of the perturbative scheme at such temperatures or the need to compute beyond next to leading 
order in ChPT, given the large relative corrections found. 

Nevertheless, the simplest realization of chiral symmetry is still provided by the much studied linear 
sigma model which possesses the convenient feature of being a renormalizable field theory, both at zero || 
and (consequently) at finite temperature Q]. A further motivation to study this model stems from recent 
theoretical results || and analyses of data J9| that seem to confirm, though not without controversy (Hi, 
that a broad scalar resonance, with a mass in the vicinity of 600 MeV -that can be identified with the 
a- meson- indeed exists. 

In addition to the above mentioned characteristics, the linear sigma model also reproduces the leading 
order modification to the pion mass in a thermal pion medium, when use is made of a systematic expansion 
in the ratio to 2 /to 2 at zeroth order, where m a , are the vacuum sigma and pion masses, respectively, 
and when treating the momentum and temperature as small quantities of order m^. In this regime, the 
effective expansion parameter becomes m 2 /47r 2 / 2 , where f v is the pion decay constant. The existence of 
such parameter allows for a controlled loop expansion of the pion self-energy from where we can extract, for 
example, the modification to the pion mass at next to leading order. 

In this work, we use the linear sigma model to compute the two-loop order modification to the pion 
propagator in a pion medium. To this end, we construct one-loop effective vertices and propagators satisfying 
the chiral Ward identities. The net result is a next to leading order calculation in the parameter to 2 /47t 2 / 2 , 
and to zeroth order in the parameter to 2 /to 2 . From the effective pion propagator, we explore the behavior 
of the pion dispersion curve at small momentum and for temperatures on the order of the pion mass. We 
also study pion damping, both in the inelastic and elastic channels. 

The work is organized as follows: In section 0, we construct an effective one-loop sigma propagator and 
one-sigma two-pion and four-pion vertices, working at zeroth order in the parameter to 2 /to 2 and satisfying 
the chiral Ward identities. In section III, we use these effective vertices and propagator to compute the 
pion self-energy at two-loop order where the effective expansion parameter turns out to be to 2 /47t 2 / 2 . We 
find the pion dispersion curve and compute also the thermal modification of the pion mass at this order. 
In section [V, we look at the damping of pions traveling in the pion medium. The main contribution to 
the damping mechanism comes from the formation of sigma resonances. We compute the mean free path 
and mean collision time for a medium pion. We also look at the damping mechanism originated in the pion 
elastic scattering processes. We finally summarize our results and conclude in section [y| A short appendix 
deals with the renormalization of the sigma model at finite temperature at two-loop order. 



II. EFFECTIVE ONE-LOOP VERTICES AND CHIRAL WARD IDENTITIES 

The Lagrangian for the linear sigma model, including only the meson degrees of freedom and after the 
explicit inclusion of the chiral symmetry breaking term, can be written as Q 

C = \ [(dnf + (daf - to 2 tt 2 - to 2 cj 2 J - A 2 /.a(a 2 + tt 2 ) - ^(a 2 + tt 2 ) 2 , (1) 

where tt and a are the pion and sigma fields, respectively, and the coupling A 2 is given by 

2 _ to 2 - to 2 
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From the above Lagrangian one obtains the Green's functions and the Feynman rules to be used in 
perturbative calculations, in the usual manner. In particular, the bare pion and sigma propagators A Tr (P), 
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Roman letters are used to denote four momenta) 



A a (Q) and the bare one-sigma two-pion and four-pion vertices r^ 2 , T^ ' are given by (hereafter, capital 



iMP)a* = 2 Sij 

iA a (Q) 



Q 2 - ml 

iT 1 ^ 1 = -2iX 2 (5 ij S kl + S ik 5 jl + 5 a 6 jk ) . (3) 

These Green's functions are sufficient to obtain the modification to the pion propagator, both at zero and 
finite temperature, at any given perturbative order. 

Alternatively, it is also possible to exploit the relations that chiral symmetry imposes among different 
n-point Green's functions. These relations, better known as chiral Ward identities (ChWI), are a direct 
consequence of the fact that the divergence of the axial current may be used as an interpolating field for 
the pion |(|. Thus, one could construct the modification to one of the above Green's functions at a given 
perturbative order and from there, build up the induced modification to other Green's functions related to 
the former by a ChWI. For example, two of the ChWI satisfied -order by order in perturbation theory- by 
the functions A 7r (P), A a (Q), r^ 2 and are 

/^jf (; 0, Pi,P 2 , P 3 ) = rg(Pi; P 2 ,P 3 )<F + r« (P 2 ;P 3 , Pi)S ik + T{ k 2 (P 3 ; P ll P 2 )5 u 

UT%(Q;0,P) = [A-\Q) A" X (P)] & , (4) 

where momentum conservation at the vertices is implied, that is Pi + P 2 + P3 =0 and Q + P = 0. The 
functional dependence of the vertices in Eqs. (Q) is such that the variables before and after the semicolon 
refer to the four-momenta of the sigma and pion fields, respectively ||. 

At one loop and after renormalization, we recall that the sigma propagator is modified by finite terms. At 
zero temperature this modification is purely imaginary and its physical origin is that a sigma particle, with a 
mass larger than twice the mass of the pion, is unstable and has a (large) non- vanishing width coming from 
its decay channel into two pions. At finite temperature the modification results in real and imaginary parts. 
The real part modifies the sigma dispersion curve whereas the imaginary part represents a temperature 
dependent contribution to the sigma width. 

The one-loop diagrams contributing to the sigma self-energy are depicted in Fig. 1 . Their explicit expres- 
sions are 

n^(g) = 6A 4 /^(Q;m 7r ,m T ), 

Ill(Q) - 18X 4 f 2 l(Q;m a , mrT ) , 
Ii% =3A\7(;m w ), 

K = -18^-/^7(; m ff ), 
mi 

Ill = -18^fij(;m a ), (5) 
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where in the imaginary-time formalism of thermal field theory (TFT), the functions X and J are defined by 

^ f ft 1 1 

X{Q;m i ,m j ) = T^J {2n)3 R2 + m 2 {K _ Q)2 + m 2 > 

r^sr^ f d 3 k 1 . , 

j(-^)=tj: j^^^ (6) 

and depend parametrically on the mass = m n , m a . Here Q — q), K — (w„,k) with uj — 2nnrT and 
u> n = 2nirT (m, n integers) being discrete boson frequencies, T is the temperature and q = |q|, k = |k|. 




w (•> (/> 



FIG. 1. Diagrams contributing to the sigma self-energy at one loop. 



Notice that the functions II* to Il£ do not develop imaginary parts. For II* the several processes described 
by a possible imaginary part are kinematically forbidden whereas for II£. to Il£ the absence of imaginary 
parts can be seen by invoking Cutkosky's rules. Moreover, as it will shortly become clear in Section III, for 
the real part of X we require only to consider its value for Q = after the analytical continuation 



iuj — > q + ie{q )q 



(7) 



which yields the time-ordered version X 1 of the function X and where e(qo) is the sign function. In this limit 
we have 



ReJ'(0; m^rrii) = -— ^5- 



dx 



whereas J is given by 



VL + x 2 

where the function / is the Bose-Einstein distribution 

/(*) 



1 + 2/ (rrnVl + x 2 ) 



+ 2/ (miVl + a; 2 ) 



(8) 



(9) 



(10) 



e x/T _ 1 ■ 

The integrals in Eqs. (||) and contain the usual T = infinities which are taken care of by the introduction 
of suitable counterterms through the renormalization procedure. The temperature-dependent terms are 
expressed as dimensionless functions of the ratio rrn/T. These are explicitly 

1 f°° dx / \ 

ReZ^Ojmi.mi) -» j ^— g / (rmy/l + x*) 
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-^h( mi /T). 



ni \ m i f dxx rf A 
27r Jo Vl + x 2 \ 



= ^(m^T), (11) 

27T Z 

where the arrows indicate only the temperature dependence of the expressions. Notice that for T ~ m^, 
g(m n /T), /i(m T /T) ~ 1. However, for to^ = m„ ~ 600 MeV, <? and /i are exponentially small and thus, in 
the kinematical regime of interest, we neglect these contributions. The remaining terms in Eq. (||) are just 
n" and H c a . Adding these up, the real part of the sigma self-energy at Q = (go, q) = becomes 

Rcn CT (0) = --^h. lh(m„/T) - j^g(mjT)i . (12) 

At this point, we make a further simplification. We notice that the ratio m 2 /A 2 / 2 can be expanded as 

#^4(1 + 4 + ...), (13) 

and that the leading term in the above expansion is of order m 2 /m 2 . We will ignore corrections of this 



order, considering the sigma particle as very heavy compared to the pion. In this approximation, Eq. (12) 
becomes 

Ren CT (0) - 6A 4 / 2 ReX t (0;TO w ,m 7r ) 
fiA 4 f 2 

= -^h(m„/T). (14) 

On the other hand, the imaginary part of II^ is expressed in terms of the imaginary part of the function I*, 
where this is given by 

ImJ*(Q; TOtt, m-n) = — > Qo + K, Q) - —> q - ie, q)} 

1 f IT / 1 — p-^+ilo, q)/T \ ~| 

= -i^ |^ 2 ) + - * i-U-tow ) } e ^ 2 - 4 ™ 2 ) > ^ 

where Q 2 = g 2 — g 2 and O is the step function and the functions a and u>± are 



/ 4m 2 



2 



IgoJ ±a(Q 2 )g . , 

w±(g ,?) = s • ( 16 ) 




FIG. 2. Diagrams contributing to the one-sigma two-pion vertex at one loop. 
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Therefore, the one-loop effective sigma propagator becomes 

1 

We now turn to the effective one-sigma two-pion vertex. The diagrams contributing at one loop order are 
shown in Fig. 2, their explicit expressions are 

SV a 12 ij (Q) = lOAV^'^Qjm^nO, 
8T\; j {Q,P) = 8X 6 f^lC(Q,P;m n ,m a ), 
SV c 12 ij (Q, P) = 24X 6 f^lC(Q, P; m a , m„) , 

ST d J j (Q) = 4\\US^I(Q; m^m a ) , 

WH j {Q) = QX 4 US^I(Q;m a ,m a ) , (18) 
where the function I is defined in the first of Eqs. ([|) and the function K, is defined by 

£(Q, P; m*, m-) = J (27r)3 K + m 2 {K _ Q? + m 2 {K _ q _ P) 2 + m 2 ■ ( 19 ) 

The expression for 5T e 12 i:> , where the internal lines are sigma propagators, is exponentially suppressed in 
the kinematical regime considered and thus we drop it from the beginning. For the functions to 
Srf^, notice also that in the kinematical regime where ra a 3> m^, T, we can pinch the internal sigma 
lines when considering the temperature dependent terms, given that the internal momentum is cut off by the 
temperature. In this approximation and after a suitable shift of the variable of integration and the analytical 
continuation in Eq. (^), it is easy to see that 

/C(Q, P; m,r, m CT ) — ► — l —I t {Q- 1 m^,m^) , 
mi 

fC(Q, P; m CT , m f ) — > — T J(; m^) , 
m% 

I(Q;m T ,m ff ) — ► -J^m-n). (20) 

Therefore, adding up the expressions in Eq. (|l8| ) after the above approximation, we have 

5T%{Q) = 2A 4 /^ j {( 5 - 4 ^) Xt (Q;m*,m*)-^ (l + J^m,)} . (21) 

As before, for the expression inside the curly brackets in Eq. (pl|), we work to zeroth order in m\jm 2 a , thus, 
the effective one-sigma two-pion vertex to one-loop order can be written as 

iT\i\Q;P 1 ,P 2 )=iT% + i8Y%{Q) 

= -2i\ 2 Ud^ [1 - SA^Qjm^m,)] , (22) 

using ir^ 2 as given in Eq. (^|). We observe that pinching the internal sigma lines and working at zeroth order 
in m^/m 2 is effectively equivalent to consider that the contributing diagrams are just those which can be 
made topologically equivalent to the diagram represented by the function m„, m,). 

Finally, we consider the effective four-pion vertex. The diagrams contributing at one loop are depicted 
in Fig. 3. Pinching the internal sigma lines and working at zeroth order in i m\jvr? a , as before, one can 
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check that the contributing diagrams are those coming from Figs. 3a to 3c, explicitly (hereafter, we omit the 
parametric dependence of the function X* on the mass rrii for the sake of a clearer notation, except when 
needed) 




FIG. 3. Diagrams contributing to the four-pion vertex at one loop 

5T^ ijkl (; P x , P 2) P 3 , P 4 ) - 2A 4 { [72*(Pj + P 2 ) + 21* (P, + P 3 ) + 2X t (P 1 + P 4 )] S ij S kl 
+ [7X t (P 1 + P a ) + 21 t (P 1 + P^ + 2I t (P 1 + P 2 )] 6* k P l 
+ [72* (Pi + Pi) + 2I t {P 1 + P 2 ) + 22*(P X + P 3 )] 5 u S lk } , 

5T b ^ kl {- P, P 2) P 3 , P 4 ) = 4A 4 { [2*(Px + P 3 ) + 2*(Pi + P 4 )] <W , 
+ [2*(Pi + P 2 ) +2*(P X + Pi)] S lk P l 
+ [l t {P l + P 2 ) +2*(Pa + P 3 )] 5"^} 

,5r^(;P 1 ,P 2 ,P 3 ,P 4 ) = -8A 4 [2*(P! +P 2 ) +X t (P 1 + P 3 ) +J t (P 1 +P i )} 

{S ij S kl +5 il 5 jk +S ik S jl } . (23) 

Adding up the expressions in Eq. ( |23[ ) and using Eq. (|^), the effective four-pion vertex up to one loop can 
be written as 

*r*/-' H (; p u p 2 , p 3 , p 4 ) = irif + isvifi; p, P 2 , P 3 , P 4 ) 

- 2iA 2 { [1 - 3A 2 2*(Pi + P 2 )] 5 i] 5 kl 
+ [l -3A 2 2*(Pi +P 3 )] 8 ik V l 

+ [1 - 3A 2 2*(Pi + P 4 )] 5 il S jk } . (24) 

It is easy to check that Eqs. (|l7|), ( p2| ) and (^) satisfy the Ward identities in Eq. (|^), this ensures that the 
approximation scheme adopted respects chiral symmetry. 

III. TWO-LOOP PION DISPERSION RELATION 

We now use the above effective vertices and propagator to construct the two-loop modification to the pion 
self-energy. The corresponding diagrams are shown in Fig. 4. Notice that diagrams Ad and 4e require knowl- 
edge of the effective, one-loop, two-sigma two-pion and three-sigma vertices, which we have not computed in 
this work. However, we can check that the leading contributions arise from diagrams 4a to 4c. To see this, 
notice that the explicit expressions for diagrams Ad and 4e are equivalent to the expressions for diagrams 
4a and 46, respectively, upon the replacement of the pion loop by a sigma loop. This results in expressions 
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exponentially suppressed and thus, in our approximation, we should ignore them. The contributing terms 
are 




O 



FIG. 4. Diagrams contributing to the pion self-energy in the effective expansion. The heavy dots denote the 
effective vertices and propagators. 



n 2 °(P) & = X 2 S ij J ±£ R2 1 ^ {5 - 91* (0) - 6A 2 Z'(P + K)} 

{6\ 2 f%\ f d 4 K 1 [ [1-3A 2 X«(0)] 2 
p (2tt)4 K* - ml \ [i_^ Xt(0 )] 



U b 2 (P) & = -\ 2 8 lj (^f-) [ 
V m a J J a 



U2iP)s -- xs I — ' l {<2.y ic- m %\ {l - ^ifi nP +K) _ a > ' 



where the subindex /? means that the integrals are to be computed at finite temperature. We now expand 
the denominators in the second and third of Eqs. (p5|) , keeping only the leading order contribution when 
considering m^, T and P as small compared to m a . Adding up the above three terms and to zeroth order 
in m 2 /m 2 where 

the pion self-energy can be written, in the imaginary-time formalism of TFT, as 

g|) [92*(0)+6Z*(P + Ji-)]} , (27) 

with K = (w n ,k) and P = (w,p). The pion dispersion relation is thus obtained from the solution to 

P 2 +m 2 +Ren 2 (P) =0, (28) 

after renormalization and the analytical continuation iuj — > p Q + it. As anticipated, to compute ReIl2(P) we 
require knowledge of the real part of X* only at Q — together with the full momentum dependence of the 
imaginary part of X* (see Eq. ( |34| ) below) . 

Notice that Eq. ( p7j ) contains temperature-dependent infinities coming from the terms involving the func- 
tion X', as well as vacuum infinities. This is an usual feature of two- loop calculations at finite temperature 
where one always encounters temperature-dependent infinities in integrals involving only the bare terms of 
the original Lagrangian. However, as it turns out, these infinities are exactly canceled by the contribution 
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from the integrals computed by using the counterterms that are necessary to introduce at one loop to carry 
the (vacuum) renormalization. The above was explicitly shown for the case of the self-energy in the 4 
theory by Kislinger and Morley JO] and for the sigma vacuum expectation value in the linear sigma model 
by Mohan |?J . We make a brief sketch of the renormalization procedure for the pion self-energy at two loops 
in the appendix to show how this cancellation occurs and refer the reader to the cited works for details. 
From now on, we concentrate on the temperature-dependent terms. 

Let us first look at the dispersion relation at leading order. After analytical continuation all the terms are 
real. The integrals involved are 

r d 3 k i i r°c dk k 2 

(27r)3#a + m2 2** J Q Ek k) 

- 9(mjT) 

rE/lM^-*®'^' (29) 

where g is the dimensionless function defined in the second of Eqs. ( |IT|) and the arrows indicate only the 
temperature dependence of the expressions. Thus, the dispersion relation results from 

[1 + 2 £ g(m w /T)] (p 2 - p 2 ) - [1 + 3 £ g(mjT)] m£ = , (30) 

with £ = m 2 /4ir 2 f 2 <C 1. For T ~ m^, g{m. K /T) ~ 1, therefore, at leading order and in the kinematical 
regime that we are considering, Eq. (Efl) can be written as 



p 2 =p 2 + ml [1 + $ g(mjT)} , (31) 

which coincides with the result obtained from ChPT ||. 

We now look at the next to leading order terms in Eq. ( |27| ) . The first of these is purely real and represents 
a constant, second order shift to the pion mass squared 

- K§)° T £/l7^^ 29Wr) ^ /r)ra "' (32) 



where h is the dimensionless function defined in the first of Eqs. (|1 1|) . The remaining term in Eq. (|27|) shows 
a non-trivial dependence on P. It involves the function S defined by 

„,„s f d 3 k 1HP + K) , s 

The sum is performed by resorting to the spectral representation of X* and (K 2 + m 2 )^ 1 . Thus, the real 
part of the retarded version of S, after analytical continuation is 

ReS r (p ,p) = - [S(iu ^p + ie,p)+S(iuj -> p - ie,p)] 



= - V 



f d 3 k [°° dk Q [°° dk' Q f, t ^,f,UV 



p -k - k' 

where V stands for the principal part of the integral. 



2tt e(^) S[k' 2 - (k - p) 2 - ml] 2Jmg(fco 1 fc) 
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Including all the terms, the dispersion relation up to next to leading order in the parameter £, for T ~ 
and in the small momentum region is obtained as the solution to 



p J =p l + { 1 + Zg(m n /T) + ^g{m v /T) [9h(mjT) - 4g(m^/T)} \ mi + £ J S(p ,p) 



where S is defined by 



5(p ,p) = -(24^ 4 )Re5 I '(po,p) 



(35) 



(36) 




FIG. 5. Pion dispersion relation obtained as the solution to Eq. ( p5[ ) for T = (upper curve). Shown is also 
the dispersion relation obtained by ignoring the term S(po,p) (lower curve). 
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FIG. 6. Temperature dependence of m w to second order in the parameter £. For comparison, the vacuum value 
is shown as a horizontal line. 



Figure 5 shows the dispersion relation obtained from Eq. ( |35[ ) for T — where we also display the 
solution without the term £ 2 S(po,p). Inclusion of this last term does not alter the shape of the curve in this 
kinematical regime. Since S(po,p) is non- vanishing for p = 0, it also contributes to the perturbative increase 
of the pion mass. Figure 6 shows the temperature dependence of m, obtained as the solution to Eq. j35|) in 
the limit when p goes to zero. Notice that the second order correction in the parameter £ has the same sign 
as the first order correction. This result is opposite to the temperature behavior of the pion mass found in 
Refs. ifl. 
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IV. PION DAMPING 



Damping of pions traveling in the pion medium is due to two possible kinds of processes: formation of 
resonances and pion elastic scattering. At temperatures on the order of the pion mass, pion elastic scattering 
is subdominant since, given their Goldston-boson nature, scattering among pions nearly disappears as the 
relative momentum approaches zero. In the linear sigma model, this feature is explicitly realized, for the 
kinematical regime discussed in this work, by the description of pion elastic scattering as a second order 
process in the small parameter £. Scattering of a pion off a given phase space element as a consequence of an 
elastic process is described by means of the imaginary part of II2 given in Eq. (^) . This can be understood 
by observing that working in the approximation where m a 3> m T , T the contributing diagrams to II2 can 
be thought of as made exclusively of pion lines, with effective one-loop vertices that involve also only pion 
lines. Thus, the rate obtained by cutting such diagrams describes a processes in which pions scatter from 
each other and keep being pions. We will discuss elastic pion scattering later on in this section. First, we 
look at resonance formation. 



FIG. 7. Relevant diagram to compute pion damping by the formation of sigma resonances. 

In the linear sigma model with only meson degrees of freedom, pions can disappear from the medium by 
forming a sigma resonance. This process is the inverse of the one describing sigma decay (except that the 
former does not happen at T=0) and thus both rates are intimately connected. The relevant diagram to 
compute resonance formation is depicted in Fig. 7 and its explicit expression is 

Hi(P) = 4\*fil(P;m r ,m a ), (37) 



where I is the function defined in Eq. (Q) that we rewrite here in terms of the quantities E* = \Jk 2 + ra 2 
and E a = \/(k — p) 2 + m 2 as 

I(P; m^)=Tj2j-^ -^L- {LUn _^ )2 + E , • (38) 



The damping rate is obtained in terms of the imaginary part of retarded version, X r of Z(P; m w , m a ), given 

by 

r/ :; /,- 

{[1 + f(\E n \) + f(\E a \)][5( Po -E n - E a ) - 5(po + E„ + E a )\ 
[/(Kl) - f(\E a \)][6(p -E n + Efj ) - 6( PQ +E W - E a )]} , (39) 



/d 3 k 1 
J^f^E~l 



after the analytical continuation iu — » p a + it. The function / above is the Bose-Einstein distribution and 



its arguments here are the absolute values of the corresponding particle energies. The integration in Eq. (39) 
can be performed analytically and the result is 
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T 
16np 



In 



1 



-E+(p ,p)/T 



1 - e -E-(p ,p)/T 



- In 



1 _ e -(E+(p ,p)+po)/T 
1 - e -(E-(p ,p)+ P0 )/T 



where the functions E± are defined by 



E±(p ,p) = b 



with 



a'(P 2 ) 



p ± a'{P 2 )p 



4m 2 
PP 1 



(40) 



(41) 



b = 



The decay rate is given in terms of Imlli by 



P Po/T 



ImIIi(po,p) 



( e Po/T _ X ) 



Po 



Pi I 



(42) 



(43) 



ImTIi contains both creation and decay rates and the factor e p °/ T /(e p °/ T — 1) in Eq. ( p3|) eliminates the 
piece describing pion creation. 

We can now compute the mean free path A for a pion traveling in the medium before forming a sigma 
resonance. This is given in terms of T> by || 



A = 



( e Po/T _ !) 



e po/T Im]li(po,p) 



(44) 



where u = p/po is the magnitude of the pion's velocity. Figure 8 shows plots of A for three different 
temperatures and a value of m a = 600 MeV as a function of the pion momentum. A reaches a maximum 
for p ~ 0.4m7r. Notice that the position of the maximum is approximately independent of T . Except for the 
location of the maximum, the behavior of the curves is both in qualitative and quantitative agreement with 
the corresponding result obtained in ChPT (see Fig. 2 Ref. @). 




FIG. 8. Mean free path as a function of momentum for a pion to travel the medium before forming a sigma 
resonance for m CT = 600 MeV. 
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FIG. 9. Mean collision time as a function of temperature for a pion to travel the medium before forming a sigma 
resonance. Shown is also the function 15/*/T J . 



We can also compute the mean collision time r between pions before forming a sigma resonance. This 
quantity is defined as the inverse of the average decay rate 
-l 



(I? 



d 3 p 
(2 



(45) 



where / is the Bose-Einstein distribution. Figure 9 shows the behavior of r as a function of the temperature 
for ra a = 600 MeV. Shown is also the quantity 15/^/T 5 . We see that r does not behave exactly as the inverse 
fifth power of the temperature p5| , though the curves are pretty close to each other in the temperature range 
considered. The pion thermal width has also been studied in Refs. |l6| p7| . 

Finally, let us discuss pion elastic scattering. The total rate, from where the cross section for elastic pion 
scattering can be obtained, is found in terms of the imaginary part of the two-loop pion self-energy II2 given 
in Eq. @ 

Imn 2 (p 0) p) = -247r 4 £ 2 Im«S> 0) p) , (46) 
where ImS r is the imaginary part of the retarded version of the function S defined in Eq. (p3), explicity J13| 



lmS r (p a ,p) = i [S(ilu — > po + ie,p) - S(iw 



IT ^( 



,Po/T 



d 3 k 
(2^)3 



dk{) 



f(ko)f(k' )S( Po -ko- k' ) 



00 dk' 



2tt 



2tt£(^) <5[fe^ — (k - 



nl] 2ImT t (fc 0l fc). 



(47) 




p/rrtir 



FIG. 10. Total pion elastic scattering rate F2 for T = m-n- as a function of momentum. Notice that this rate is 
about two orders of magnitude smaller than Fi, given that it is proportional to £ 2 . 
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FIG. 11. One and two-loop diagrams and counterterms for the computation of the pion-self energy described by 
the effective Lagrangian in Eq. (|5ll). The temperature-dependent infinities are exactly canceled. 

The total reaction rate is given by 

T2(Po,p) = - — ImIT 2 (p ,p) , (48) 
Po 

where po should be taken as the solution to Eq. (pl|). The integration in Eq. can only be performed 
numerically Figure 10 shows a plot of F 2 as a function of p for T = m^. As anticipated, this rate is about 
two orders of magnitude smaller than T±, given that it is proportional to £ 2 . 



V. SUMMARY AND CONCLUSIONS 



In conclusion, working in the linear sigma model at finite temperature, we have found effective one- 
loop sigma propagator and one-sigma two-pion and four-pion vertices satisfying the ChWI and therefore 
respecting chiral symmetry, when maintaining only the zeroth order terms in an expansion in the parameter 
to 2 /to 2 . We have used these objects to compute the two- loop order correction to the pion propagator in 
a pion medium for small momentum and for T ~ and showed that the linear sigma model yields the 
same result as ChPT at leading order in the parameter £ = to 2 /47t 2 / 2 . This result was to be expected since 
the kinematical regime we consider is that where the temperature, the pion momentum and the pion mass 
are treated as small quantities, such as in the case of ChPT. The existence of the parameter £ allows for a 
controlled perturbative calculation of physical quantities in the above mentioned kinematical regime, such 
as the thermal modification of the pion mass. We have shown that, contrary to the result in Rcfs. the 
two-loop order correction to the pion mass is proportional to £ 2 and that this correction is of the same sign 
as the one-loop correction. The shape of the dispersion curve is not significantly altered in the kinematical 
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regime considered. We should also note that our result is more general than the one obtained in Ref. JL8), 
where the computation of the two-loop pion dispersion relation is carried from the onset in the weak coupling 
limit, that is to say, for A 2 <C 1. In our approach, the first step can be thought of as valid for arbitrary values 
of A 2 , this is why the effective one-loop vertices and propagator found depend explicitly on A 2 . However, 
when it comes to the computation of a physical quantity such as the pion dispersion relation, we have shown 
(see Eq. J26|)) that the contribution from A 2 enters the result always in the combination 



showing that for m a > m„, the leading contribution is independent of m a . 

We have also studied pion damping both in the elastic and inelastic channels. We have computed the 
momentum dependence of the mean free path for a pion before forming a sigma resonance and found that 
the shape of the curves coincide both quantitatively and qualitatively with the result obtained from ChPT 
when using a value of m a — 600 MeV. The mean collision time of a pion before forming a sigma resonance 
as a function of T is however close but not quite proportional to the inverse fifth power of the temperature. 
It remains to include the effects of a possible high nucleonic density in the same kinematical regime. 

Regarding the applicability of the linear sigma model as a tool to study low temperature QCD, let us 
stress that at such energies, any theory whose aim is to describe such degrees of freedom is an effective one. 
The essential ingredient, common to all of them is chiral symmetry, that is to say, the consequences arising 
from the smallness of the pion mass. Therefore, as long as chiral symmetry is respected and the temperature 
is small enough (on the order of the pion mass) so as to not excite higher resonances, all such theories will 
yield equivalent results. We have shown explicitly that this is the case since, contrary to common belief, 
the one-loop computation of the pion dispersion curve yields the same result as the one obtained by using 
ChPT, provided we work in the same kinematical regime. 

We should also mention that for temperatures on the order of 150 MeV, it might be important to include 
the effects introduced by other degrees of freedom, such as the p meson. In this sense, the quantities computed 
in this work at such temperatures must be regarded as examples of the results obtained by considering only 
the lowest lying meson degrees of freedom and by no means as the final answer to the pion propagation 
properties. 

Finally, it is interesting to note that the formalism thus developed could be employed to explore the 
behavior of the pion dispersion curve in the large momentum region where, in principle, the lowest order 
ChPT Lagrangian cannot be used. These issues will be treated in a follow up work Jig]. 
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APPENDIX: TWO-LOOP RENORMALIZATION AT FINITE TEMPERATURE 



The purpose of this appendix is to show how the cancellation of the temperature-dependent infinities 
appearing in two loop computations is achieved. For this matter, we concentrate on the two-loop pion self- 
energy given in Eq. (p7|). Notice that the first and third terms can be nicely combined to give a numerical 
factor of 3. Moreover, we carry the renormalization procedure at P = where the approximation used 
renders the result exact. Thus, Eq. (27) can be written as 

iw> = (g)r£/|^_±_{ 3 _ (g) mo)+ ^ {P + K)] ) . m 

The two-loop self-energy expressed in Eq. (JSTj) can be formally obtained by means of the effective Lagrangian 

C = \{d^f-\ m l^-^(^) 2 , (51) 

where a = 6(m^./2/^). m T and are taken as the renormalized values of the pion mass and decay constant, 
respectively. The factor 6 comes from considering the interaction of like-isospin pions in the vertex 

irf = -2i (5 ij 6 kl + S ik S jl + 8 a V k ) . (52) 

At one-loop, there are two (temperature-independent) infinities coming from diagrams a and b in Fig. 11. 
Using the effective Lagrangian in Eq. ( pi] ) and dimensional regularization, the counterterms needed to cancel 
these infinities, represented by diagrams c and d in Fig. 11, are given explicitly by 

Smt = 3 



2 ft) IBttV 



2 \ 2 



fa=3x36 UfJ ife' (53) 

where e = d — 4, d being the number of dimensions, and the factor 3 in the second of Eqs. (^) takes care of 
the three crossed channels represented by Fig. 116, since the infinity in each of these is the same, /i is the 
usual mass parameter introduced in dimensional regularization. 

Going to two loops, temperature-dependent infinities appear. Diagram e in Fig. 11, which explicit ex- 
pression is given by the second term in Eq. ( p0| ) gives rise to two infinite, temperature-dependent terms, 
coming from the product of the vacuum infinity in I t (0) times the temperature-dependent term in the inte- 
gration and the vacuum infinity in the integration times the temperature-dependent piece in 2"*(0). These 
are explicitly 

a — *' (#) 2 i&@)^>. 



2 



im: 



2 



* = - 9 [W ^ (0) ' (54) 



where g(m n /T) is the dimensionless function defined in the second of Eqs. ( pT| ) and 2y.(0) represents the 
temperature-dependent part of the function X 1 (0) . 

Diagram / in Fig. 11, which explicit expression corresponds to the third term in Eq. (|50|), gives rise to the 
temperature-dependent infinity 
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y "- m iW ^k{m s{m ' /T> - (55) 

To complete the calculation, we need to compute the integrals involving the one-loop counterterms. Diagram 
g in Fig. 11 gives the explicit temperature-dependent infinity 

2 

2 " V 2 /^/ iB^e V 271-2 



% ^74(0). (57) 



whereas diagram h in Fig. 11 gives 

2 

,2/2 J 16^ 2 e 

We can now see that 3^5 in Eq. (|57|) cancels 3^2 in the second of Eqs. ( |54| ) and that 3^1; 3^3 and 3^4 in 
Eqs. (p3), ( |55| ) and (|5^), cancel among themselves. Therefore, the result at two loops is free form temperature- 
dependent infinities. Similar arguments apply to other Green's functions, see Refs. @@ for additional 
details. 
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